Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 12/14/05 



REFINED PARAMETERS OF THE PLANET ORBITING HD 189733 



o 
o 

(N 



(N 



G. A. Bakos^'^, H. Knutson\ F. Pont^, C. Moutou^, D. Charbonneau\ A. Shporer*, F. Bouchy'' "\ M. Everett'' 
C. Hergenrother' , D. W. Latham\ M. 



Mayor'", T. MAZEH^ R. W. Noyes\ 

Draft version February 5, 2008 

ABSTRACT 



D. QuELOZ"', A. Pal and S. Udry"' 



We report on the BVRI multi-band follow-up photometry of the transiting extrasolar planet 
HD 189733b. We revise the transit parameters and find planetary radius i?p = 1.154 ± 0.032i?j 
and inclination ip — 85.79 ± 0.24°. The new density Igcm"'^) is significantly higher than the 
former estimate 0.75gcm~'^); this shows that from the current sample of 9 transiting planets, only 
HD 209458 (and possibly OGLE-lOb) have anomalously large radii and low densities. We note that 
due to the proximity of the parent star, HD 189733b currently has one of the most precise radius 
determinations among extrasolar planets. We calculate new ephemerides: P = 2.218573 ± 0.000020 
days, To = 2453629.39420 ± 0.00024 (HJD), and estimate the timing offsets of the 11 distinct transits 
with respect to the predictions of a constant orbital period, which can be used to reveal the presence 
of additional planets in the system. 

Subject headings: stars: individual: HD 189733 - planetary systems 
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1. INTRODUCTION 

HD 189733 is one of nine currently known main se- 
quence stars orbited by a transiting giant planet. The 
system is of exceptional interest because it is the clos- 
est known transiting planet (D = 19.3pc), and thus is 
amenable to a ho st of follow-up obser vations. Th e dis- 
covery paper by iBouchv et al.l l)2005j) (hereafter lB05fl 
derived the key physical characteristics of the planet, 
namely its mass (1.15 ± 0.04Afj) and radius (1.26 ± 
0.03i?j), based on radial velocity observations of the 
star made with the ELODIE spectrograph at the 1.93m 
telescope at Observatoire de Haute Provence (OHP), to- 
gether with photometric measurements of one complete 
and two partial transits made with the 1.2m telescope 
also at OHP. With these parameters HD 189733 b had a 
large radius comparable to HD 209458b ( Laughli n et al.l 
|2QQ3); s-nd a density roughly equal to that of Saturn 
{p ~ 0.7gcm~^). 

Determining precise radii of extrasolar planets in ad- 
dition to their mass is an im portant foc us of exqplane t 
research (see e.g. IBouchv et al. 2004; Tor res et al.ll2004j) . 
because the mean density of the planets can shed light 
on their internal str ucture and evolution. According to 
iBaraffe et al.l (|2Qn3)> the radii of all known extrasolar 
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planets are broadly consistent with models, except for 
HD 209458b. This planet with its large radius and low 
density {p ~ 0.33gcm^'^) has attracted considerable in- 
terest, and various mechanisms involving heat deposition 
beneath the surface have been suggested (jLaughlin et alJ 
120051 and references therein). An additional motivation 
for obtaining accurate planetary radii is proper interpre- 
tation of follow-up data, notably secondary eclipse and 
reflected light observations. This is of particular rele- 
vance to HD 189733b, which h as been recently obs erved 
by the Spitzer Space Telescope l|Deming et al.l20d^ . and 
where the brightness temperature depends on the radius 
ratio of the planet to the star. 

Both by extending the current, very limited sample of 
transiting exoplanets, and by precise determination of 
the physical parameters it will become possible to refine 
theoretical models and decide which planets are "typ- 
ical". Close-by, bright stars, such as the host star of 
HD 189733 are essential in this undertaking. The OGLE 
project (Udalski ct al. 2002a b e) and fo llow-up obser- 
vations (e.g. Konacki ct al. 2004; Mout ou et all 120041: 
made a pivotal contribution to the cur- 
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rent sample by the discovery of more than half of the 
known transiting planets. Follow-up observations, how- 
ever, are cumbersome due to the faintness of the targets, 
and required the largest available telescopes. The typ- 
ical errors of mass and radius for these host stars are 
^ O.O6M0 and ~ O.15i?0, and the corresponding errors 
in planetary parameters are ^ 0.13Mj and ^ 0.12_Rj, 
respectively. However, for planets orbiting bright stars 
in the solar neighborhood, errors at the level of a few 
percent can be reached for both the mass and radius. 

In this paper we report a number of follow-up photo- 
metric measurements of HD 189733 using six telescopes 
spaced around the world. Together with the original 
OHP photometry, we use these measurements to deter- 
mine revised values for the transit parameters, and give 
new ephemerides. First we describe the follow-on pho- 
tometry in detail (32)), followed by the modeling which 
leads to the revised estimate of the planetary radius (^SJ, 
and we conclude the paper in 21 
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TABLE 1 

Summary of instruments used in the observing campaign of HD 189733. 



Site 


Longitude 


Latitude 


Alt. 


Telescope 


Diam. 


Detector 


Pxs 


Trd 


FOV 












(meters) 




(meters) 




C'/pix) 


(sec) 




OHP 


05° 30' 


E 


43° 55' 


N 


650 


OHP1.2 


1.2 


SITe IK X IK 


0.69 


90 


11.77' 


FLWO 


110° 53' 


W 


31°41' 


N 


2350 


FLWO 1.2 


1.2 


Fairchild 4K X 4K 


0.34 


12 


23' 


FLWO 


110° 53' 


W 


31°41' 


N 


2345 


HAT-5 


0.11 


Thomson 2K X 2K 


14.0 


10 


8.2° 


FLWO 


110° 53' 


W 


31°41' 


N 


2345 


HAT-6 


0.11 


Thomson 2K X 2K 


14.0 


10 


8.2° 


FLWO 


110° 53' 


W 


31°41' 


N 


2345 


TopHAT 


0.26 


Marconi 2K X 2K 


2.2 


40 


1.29° 


Mauna Kea 


155° 28' 


W 


19°49' 


N 


4163 


HAT-9 


0.11 


Thomson 2K X 2K 


14.0 


10 


8.2° 


Wise 


34°35' 


E 


30° 35' 


N 


875 


Wisel.O 


1.0 


Tektronics IK X IK 


0.7 


40 


11.88' 



2. OBSERVATIONS AND DATA REDUCTION 

We organized an extensive observing campaign with 
the goal of acquiring muhi-band photometric measure- 
ments of the transits of HD 189733 caused by the hot 
Jupiter companion. Including the discovery data of lBOfil 
that were obtained at OHP, altogether four sites with 
seven telescopes contributed data to two full and eight 
partial transits in Johnson B, V, R, I and Sloan r photo- 
metric bandpasses. The sites and telescopes employed 
are spread in geographic longitude, which facilitated 
gathering the large number (close to 3000) of individual 
data spanning 2 months. 

The following telescopes were involved in the photo- 
metric monitoring: the Im telescope at the Wise Obser- 
vatory, Israel; the 1.2m telescope at OHP; the 1.2m tele- 
scope at Fred Lawrence Whipple Observatory (FLWO) 
of the Smithsonian Astrophysical Observatory (SAO); 
the 0.11m HAT-5 and HAT-6 wide field telescopes plus 
the 0.26m TopHAT telescope also at FLWO; and the 
0.11m HAT-9 telescope at the Submillimeter Array site 
at Mauna Kea, Hawaii. An overview of the sites and 
telescopes is shown in Tabled 

A summary of the observations is shown in Table [3 
The telescopes are identified by the same names as in 
Tabled The transits have been numbered starting with 
the discovery data Ntr = 0, and are identified later in 
the text using this number. In the following subsections 
we summarize the observations and reductions that are 
specific to the sites or instruments. 

2.1. Observations by OHP 1.2m telescope 

These obs ervati ons and their reduction were already 
described in lB05i Summarizing briefly, the 1.20m f/6 
telescope was used together with a IK x IK back- 
illuminated CCD having 0.69"/pix resolution. Typical 
exposure times were 6 seconds long, followed by a 90 sec- 
ond readout. The images were slightly defocused, with 
FWHM«2.8". 

Full-transit data were obtained in Johnson B-band un- 
der photometric conditions for the Ntr = transit. This 
is shown by the "OIBEO" flag in Table |21 indicating 
that the Out-of-transit part before the Ingress, the Bot- 
tom, the Egress, and the Out-of-transit part after egress 
all have been observed. This is an important part of 
the combined dataset, as it is the only full transit seen 
in B-band. In addition, partial transit data were ob- 
tained for the Ntr = 4 event using Cousin's R-band fil- 
ter (Re) under acceptable photometric conditions, and 
for the Ntr = 5 event using the same filter under non- 
photometric conditions. The frames were subject to bias. 



dark and fiatfield calibration procedure followed by cos- 
mic ray removal. Aperture photometry was performed in 
an aperture of 9.6" using the daophot (Stetson 198^ 
package. 

The B-band light-curve published in the IB05I discov- 
ery paper used the single comparison star HD 345459. 
This light-curve suffered from a strong residual trend, 
as suggested by the ^ O.Olmag difference between the 
pre- and post-transit sections. This trend was probably 
a consequence of differential atmospheric extinction, and 
was removed by a linear airmass correction, bringing the 
two sections to the same mean value. This ad-hoc cor- 
rection, however, may have introduced an error in the 
transit depth. In this paper, we used six comparison 
stars in the field of view (selected to have comparable 
relative flux to HD 189733 before and after transit). A 
reference light-curve was built by co-adding the normal- 
ized flux of all six stars, and was subtracted from the 
normalized light-curve of HD 189733. The new reduc- 
tion shows a residual OOT slope 4.2 times smaller than 
in the earlier reduction. The resulting transit depth in 
B-band is decreased by about 20% compared to the dis- 
covery data. This illustrates the large contribution of 
photometric systematics that must be accounted for in 
this kind of measurement. The R-band data set is not 
as sensitive to the extinction effect as the B-band, hence 
the selection of comparison stars has a minimal impact 
on the shape of the transit curve. The B-band light-curve 
is shown on panel 6 of Fig.Q] the R-band light-curves are 
exhibited in Fig. [5] 

2.2. Observations by the FLWO 1.2m telescope 

We used the FLWO 1.2m telescope to observe the full 
transit of Ntr = 6 in Sloan r band. The detector was 
Keplercam , which is a single chip 4K x 4K CCD with 
15/im pixels that correspond to 0.34" on the sky. The 
entire field-of-view is 23'. The chip is read out by 4 am- 
plifiers, yielding a 12 second readout with the 2x2 bin- 
ning we used. The single-chip design, wide field-of-view, 
high sensitivity and fast readout make this instrument 
well-suited for high-quality photometry follow-up. 

The target was deliberately defocused in order to al- 
low longer exposure times without saturating the pixels, 
and to smear out the inter-pixel variations that may re- 
main after fiatfield calibration. The intrinsic FWHM was 
~ 2", which was defocused to ~ 10". While conditions 
during the transit were photometric^^, there were partial 

This was confirmed from the all-sky webcamera movies taken 
at the Multiple-Mirror Telescope (MMT) that are archived on a 
nightly basis tjhttp: / /skycam. mmto.arizona.edu/ I 
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TABLE 2 

Summary of HD 189733 observations. 



Telescope 


Filter 


Ntr 


Epoch 


Date 


Transit 


Cond. 


TOOT 


^sys 


Cad. 


A r//i 

Ap["] 










(UT) 






(mmag) 


(mmag) 


(sec) 


(") 


OHP1.2 


B 







2005-09-1 5 


GIBEQ 


5 


2.6 


1.3 


86 


10 


Wispi n 

V V A- *\j 


B 


4 


^3638 3 


2005-09-24 


-IBE- 


4 






42 


5 


OHP1.2 




4 


53638 3 


2005-09-24 


— BEQ 


4 


3.0 


1.2 


95 


10 


OHP1.2 




5 


53640 5 


200,5-09-26 


01 


3 


6.8 


2.4 


95 


10 


FTWOI 2 




g 


53642 7 


2005-09-29 


OIBEQ 


4b 


2.6 


0.5 


17 


20 


HAT-5 


Ic 


6 


53642.7 


2005-09-29 


OIBEQ 


4b 


4.4 


1.3 


135 


42 


HAT-6 


Ic 


6 


53642.7 


2005-09-29 


OIBED 


4b 


4.1 


1.2 


108 


42 


TopHAT 


V 


6 


53642.7 


2005-09-29 


OIBEQ 


4b 


4.6 


3.0 


70 


10 


HAT-9 


Ic 


7 


53644.9 


2005-10-01 


OIB— 


4 


4.6 


1.2 


99 


42 


HAT-9 


Ic 


16 


53664.9 


2005-10-21 


OIB— 


4 


4.3 




100 


42 


TopHAT 


V 


19 


53671.6 


2005-10-28 


eO 


4 


5.3 




106 


10 


HAT-5 


Ic 


19 


53671.6 


2005-10-28 


EQ 


4 


4.6 




103 


42 


HAT-5 


Ic 


20 


53673.8 


2005-10-30 


Olb— 


5 


3.3 


0.9 


85 


42 


Wisel.O 


B 


22 


53678.2 


2005-11-03 


~bEQ 


3 


5.5 


1.1 


49 


10 


TopHAT 


V 


24 


53682.6 


2005-11-08 


Olb— 


2 


5.4 


2.6 


108 


10 


HAT-9 


Ic 


29 


53693.7 


2005-11-19 


-IBEO 


5 


6.6 


1.1 


90 


42 



Note. — The table summarizes all observations that were part of the observing campaign described in this 
paper. Not all of them were used for refining the ephemerides or parameters of the transit — see Tabled and 
Tabic 1^ for reference. Ntr shows the number of transits since the discovery data. Epoch and Date show the 
approximate time of mid-transit. The Transit column describes in a terse format which parts of the transits were 
observed; Out-of- Transit (OOT) section before the transit, Ingress, Bottom, Egress and OOT after the transit. 
Missing sections are indicated by The Conditions column indicates the photometric conditions on the scale 

of 1 to 5, where 5 is absolute photometric, 4 is photometric most of the time with occasional cirrus/fog (relative 
photometric), 3 stands for broken cirrus, and 2 for poor conditions. Column (Tqot gives the rms of the OOT 
section at the Cadence shown in the next column. If the transit was full, (Joot was computed separately from 
the pre- and post-transit data, and the smaller value is shown. Column Gsys shows the estimated amplitude of 
systematics (for details, see . Ap shows the aperture used in the photometry in arcseconds. 

^Sloan r filter 

^Conditions were non-photometric before the transit (on the initial part of the OOT), then they became photo- 
metric for the entire duration of the transit, and deteriorated after the transit at the end of the OOT. 



clouds before and after. The focus setting was changed 
twice during the night; first when the clouds cleared, and 
second, when the seeing improved. In both cases the rea- 
son was to keep the signal level within the linear response 
range of the CCD. We used a large enough aperture 
that these focus changes did not affect the photometry. 
All exposures were 5 seconds in length with 12 seconds 
of readout and overhead time between exposures. We 
observed the target in a single band so as to maximize 
the cadence, and to eliminate flatfielding errors that may 
originate from the imperfect sub-pixel re-positioning of 
the filter-wheel. Auto-guiding was used to further mini- 
mize systematic errors that originate from the star drift- 
ing away on the CCD chip and falling on pixels with 
different (and not perfectly calibrated) characteristics. 

Reduction and photometry — All images were reduced in 
the same manner; applying overscan correction, subtrac- 
tion of the two-dimensional residual bias pattern and cor- 
rection for shutter effects. Finally we flattened each im- 
age using a combined and normalized set of twilight sky 
flat images. There was a drift of only ~ 3" in point- 
ing during the night, so any large-scale flatfielding errors 
were negligible. 

To produce a transit light curve, we chose one image 
as an astrometric reference and identified star centers for 
HD 189733 and 23 other bright and relatively uncrowded 
stars in the field. We measured the flux of each star 
around a fixed pixel center derived from an astrometric 
fit to the reference stars, in a 20" circular aperture using 
daophot/phot within iraf^^ (|Todvl 11980. Il993() and 

IRAF is distributed by the National Optical Astronomy Ob- 



estimated the sky using the sigma-rejected mode in an 
annulus defined around each star with inner and outer 
radii of 33" and 60" respectively. 

We calculated the extinction correction based on a 
weighted mean flux of comparison stars and applied this 
correction to each of our stars. We iteratively selected 
our comparison stars by removing any that showed un- 
usually noisy or variable trends in their differential light 
curves. Additionally, a few exposures in the beginning 
and very end of our observing sequence were removed 
because those observations were made through particu- 
larly thick clouds. The resulting light curve represents 
the observed counts for the star corrected for extinction 
using a group of 6 comparison stars within 6' separation 
from HD 189733. The light-curve is shown on panels 3 
and 4 of Fig. [l] 

2.3. Observations by HATNet 

An instrument de scription of the wide-field HAT tele- 
scopes was given in iBakos et al.l l)2002l 120041) . Here we 
briefly recall the relevant system parameters. A HAT 
instrument contains a fast focal ratio (f/1.8) 0.11m di- 
ameter Canon lens and Peltier-cooled CCD with a front- 
illuminated 2K X 2K chip having 14/Ltm pixel size. The 
resulting FOV is 8.2° with 14" pixe l scale. Using a psf- 
broadening technique ( Ba kos et alJ 12004') . careful cali- 
bration procedure, and robust differential photometry, 
the HAT telescopes can achieve 3mmag precision (rms) 
light-curves at 300s resolution for bright stars (at / « 8) . 

servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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The HAT instruments are operated in autonomous mode, 
and carry out robotic observations every clear night. 

We have set up a longitude-separated, two-site network 
of six HAT instruments, with the primary goal being 
detection of planetary transits in front of bright stars. 
The two sites are FLWO, in Arizona, the same site where 
the 1.2m telescope is located f i)2.2|) . and the roof of the 
Submillimeter Array atop Mauna Kea, Hawaii (MK). 

In addition to the wide-field HAT instruments, we de- 
veloped a dedicated photometry follow-up instrument, 
called TopHAT, which is inst alled at FLWO. A br i ef sys- 
tem description was given in iCharbonneau et al.l l|20flf]f) 
in context of the photometry follow-up of the HD 149026 
planetary transit. This telescope is 0.26m diameter, f/5 
Ritchey-Cretien design with a Baker wide-field corrector. 
The CCD is a 2K X 2K Marconi chip with 13.5^m pixel 
size. The resulting FOV is 1.3° with 2.2" pixel resolution. 
Similarly to the HATs, TopHAT is fully automated. 

Selected stations of the HAT Network, along with 
TopHAT, observed one full and six partial transits of 
HD 189733 (for details, see Table 121. Observing condi- 
tions of the full-transit event at FLWO at Nt^ = 6 have 
been summarized in t l2.2l This transit was observed by 
HAT-5 and HAT-6 (both in I-band), and by TopHAT 
(V-band). The partial transit observations at numerous 
later epochs included HAT-5 (FLWO, I-band), HAT-9 
(MK, I-band) and TopHAT (FLWO, V-band). Typical 
exposure times for the wide-field instruments were 60 to 
90 seconds with 10 second readout. TopHAT exposures 
were '^12sec long with up to 40 second readout and down- 
load time. All observations were made at slight defocus- 
ing and using the psf-broadening technique. The stellar 
profiles were 2.5pix (35") and 4.5pix (9.9") wide for the 
HATs and TopHAT, respectively. Although we have no 
auto-guiding, real-time astrometry was performed after 
the exposures, and the telescope's position was kept con- 
stant with 20" accuracy. 

Reduction and photometry — All HAT and TopHAT 
images were subject to overscan correction, two- 
dimensional residual bias pattern and dark subtraction, 
and normalization with a master sky-fiat frame. Bias, 
dark and sky-flat calibration frames were taken each 
night by each telescope, and all object frames were cor- 
rected with the master calibration images that belonged 
to the specific observing session. Saturated pixels were 
masked before the calibration procedure. 

We used the 2MASS All-Sky Catalog of Point Sources 
ijSkrutskie et al.ll2000t ICutri et alJl2003|) as an input as- 
trometric catalogue, where the quoted precision is 120 
mas for bright sources. A 4th order polynomial fit was 
used to transform the 2MASS positions to the reference 
frame of the individual images. Typical rms of the trans- 
formations was 700 mas for the wide- field instruments, 
and 150 mas for TopHAT. 

Fixed center aperture photometry was performed for 
all these stars. For the wide-field HAT telescopes we 
used an Tap = 3 pixel (42") aperture, surrounded by an 
annulus with inner and outer radii of ri = 5 pix (70") 
and r2 = 13 pix (3'), respectively. For TopHAT, the best 
aperture was rap = 5 pix (10.8") with ri = 13 pix (29") 
and r2 = 21 pix (46"). The apertures were small enough 
to exclude any bright neighboring star. 

A high quality reference frame was selected for the 



wide-field HAT telescopes from the Mauna Kea HAT- 
9 data, and separately for TopHAT. Because the HAT 
wide field instruments are almost identical, we were able 
to use the HAT-9 reference frame to transform the in- 
strumental magnitudes of HAT-5 and HAT-6 data to a 
common system. For this, we used 4th order polynomi- 
als of the magnitude differences as a function of X and 
Y pixel positions. In effect, we thereby used ~ 3000 
and ~ 800 selected non-variable comparison stars for the 
HATs, and TopHAT, respectively. This contributes to 
the achieved precision, which is only slightly inferior to 
the precision achieved by the bigger diameter telescopes. 

The amount of magnitude correction for HD 189733 
between the reference and the individual images is shown 
in the l^Mf-xt column of Table |3| The same table also 
indicates the rms of these magnitude fits in the Omfit 
column. Both quantities are useful for further cleaning 
of the data. Because HD 189733 is a bright source, it was 
saturated on a small fraction of the frames. Saturated 
data-points were fiagged in the light-curves, and also de- 
selected from the subsequent analysis (flagged as "C" in 
Table EJ. 

After cleaning outliers by automatically de-selecting 
points where the rms of the magnitude transformations 
was above a critical threshold (typically 25mmag) the 
light-curves reached a precision of ~4mmag at 90 second 
resolution for both the HATs and TopHAT. Full-transit 
data are shown in panels 1, 2 and 5 of Fig.^ and partial- 
transit data are shown in Fig. El 

2.4. Observations by the Wise 1.0m telescope 

The Wise Im f/7 telescope was used to observe the 
Ntr = 4 and Ntr = 22 transits in B-band. The CCD 
was a IK x IK Tektronics chip with 24/im pixel size that 
corresponds to 0.696"/pix resolution on the sky, and a 
FOV of 11.88'. The photometric conditions were accept- 
able on both nights, with FWHM«2". Auto-guiding was 
used during the observations. Frames were calibrated in 
a similar manner to the FLWO 1.2m observations, using 
twilight flats, and aperture photometry was performed 
with DAOPHOT. 

Unfortunately, out-of-transit (OOT) data of the first 
transit {Ntr = 4) (which was also observed from OHP 
in R-band) are missing, so it is impossible to obtain use- 
ful normalization or to apply extinction correction to the 
transit curve. The second transit {Ntr = 22) was pro- 
cessed using an aperture of 10" encircled by an annulus 
with inner and outer radii of 15" and 25", respectively. 
Seven comparison stars were used, all of them bright, iso- 
lated and far from the boundary of the FOV. Extinction 
correction, derived from the OOT points only, was ap- 
plied to the resulting stellar light-curve. The final curve 
of this transit is plotted in Fig. |21 

2.5. The resulting light-curve 

All photometry originating from the individual tele- 
scopes which contains significant OOT data has been 
merged, and is presented in Table |3 We give both the 
ratio of the observed flux to the OOT flux of HD 189733 
( "FR" ) , and magnitudes that are very close to the stan- 
dard Johnson/Cousins system ("Mag"). Due to the dif- 
ferent observing conditions, instruments, photometry pa- 
rameters (primarily the aperture) and various systematic 
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effects (changing FWHM), the zero-point of the obser- 
vations were slightly offset. Even for the same instru- 
ment, filter-setup, and magnitude reference frame, the 
zero-points in the flat OOT section were seen to dif- 
fer by 0.03mag. The offset can be explained by long- 
term systematic variations and by intrinsic variation of 
HD 189733. 

In order to correct for the offsets, for each transit ob- 
servation (as indicated by Ntr in Table we calculated 
both the median value from the OOT section by reject- 
ing outliers, and also the rms around the median. The 
OOT median was used for two purposes. First, we nor- 
malized the flux values of the given light-curve segment 
at Ntr , which are shown in the "FR" (flux-ratio) column 
of Table 121 Second, we shifted the magnitudes to the 
standard system in order to present reasonable values in 
the "Mag" column. For the standard system we used the 
Hipparcos values, except for R-band, which was derived 
by assuming R~I ^ 0.48 from Cox (2000). 

The formal magnitude errors that are given in the 
"Merr" column are based on the ph oton-noise of th e 
source and the background- noise (e.g. lNewberrvlll99H) . 
They are in a self-consistent system, but they un- 
derestimate the real errors, which have contri butions 
from other noise factors, such as i) scintillation ijYound 
Il^ jl7: Gilli l and ^ Brown 1988), ii) calibration frames 
ijNewberrvi il99ll) . iii) magnitude transformations de- 
pending on the reference stars and imperfectly corrected 
extinction (indicators of this error source are the AM^xt 
extinction and the amfu rnis of extinction corrections in 
Table O. Because it is rather difficult to calculate these 
factors, we assumed that the observed rms in the OOT 
section of the light-curves is a relevant measure of the 
overall noise, and used this to normalize the error esti- 
mates of the individual flux-ratios (column "FRerr" , see 
later 

2.6. Merger analysis 

HD 189733 has a number of faint, close-by neighbors 
that can distort the light-curve, and may bias the derived 
physical parameters. These blends can have the following 
second-order effects: i) the measured transit will appear 
shallower, as if the planetary radius was smaller, ii) the 
depth and shape of the transit will be color-dependent 
in a different way than one would expect from limb- 
darkening models, iii) differential extinction can yield an 
asymmetric light-curve, iv) variability of a faint blend 
can influence the observed light-curve. Our goal was to 
calculate the additional flux in the various apertures and 
bandpasses shown in Table [51 and correct our observed 
flux-ratios (Table 13 column FR) to a realistic flux-ratio 
(FR 

The 2MASS point-source catalogue ijCutri et alJl2003(l 
lists some 30 stars within 45", which is the aperture 
used at the HAT-5,6,9 telescopes, 5 stars within 20" 
(FLW01.2m), and 3 stars within 15", which may af- 
fect the measurements of the 10" apertures of 0HP1.2, 
Wisel.O and TopHAT (apertures are hsted in Table EJ. 
To check the reality of listed blends, and to search for ad- 
ditional merger stars, we inspected the following sources: 
the Palomar Observatory Sky Survey (POSS I) red plates 
(epoch 1951), the Palomar Quick-V Survey (QuickV, 
epoch 1982), the Second Palomar Sky Survey (POSS 
II) plates (epoch 1990 - 1996), the 2MASS J, H and 



Ks scans (epoch 2000), and our own images. We can 
use the fact that HD 189733 is a high proper motion 
star with velocity of ~ 0.25" /yr pointing South. It was 
~ 13" to the North on the POSS-I plates, and 4" N on 
POSS-II, thus we can check its present place when it was 
not hidden by the glare of HD 189733. The analysis is 
complicated by the saturation, diffraction spikes, and the 
hmited scan resolution (1.7"/pix) of POSS-I, but we can 
confirm that there is no significant source at the epoch 
2005 position of HD 189733 down to '--^4mag fainter in R- 
band. The reality of all the 2MASS entries was double 
checked on the POSS frames. 

There are only two additional faint sources that are 
missing from the 2MASS point-source catalogue, but 
detected by our star-extraction on the 2MASS J, H 
and K scans; the first one at a = 20'*00"45.12^ S = 
-h22°42'36.5" and the second at a = 20'*00"43.20^ 
6 = -|-22°42'42.5". We made sure these sources are not 
filter-glints or persistence effects on the 2MASS scans; 
they are also visible on the POSS frames. Their instru- 
mental magnitude was transformed to the J,H,_ftr5 system 
using the other stars in the field that are identified in the 
point source catalogue. 

A rough linear transformation was derived between 
the 2MASS J, H and Ks colors and Johnson/Cousins 
B, V, R and I by cross-identification of ^ 450 Landol3 
l)199^ standard stars, and performing linear regression. 
The uncertainty in the transformation can be as large as 
O.lmag, but this is adequate for the purpose of estimat- 
ing the extra flux (in BVRI), which is only about a few 
percent that of HD 189733. 

We flnd that the extra flux in a 45" aperture is 
S = 1.012, 1.016, 1.018 and 1.022 times the flux of 
HD 189733 in B, V, R and I-bands, respectively. The 
dominant contribution comes from the red star 2MASS 
20004297-1-2242342 at 11.5" distance , which is^4.5m ag 
fainter. This star has been found l(Bakos et alJl200^ to 
be a physical companion to HD 189733 and thus may 
also be called HD 189733B (not to be confused with 
HD 189733b). For the 10" aperture we assumed that half 
the flux of HD 189733B is within the aperture. The same 
S flux contribution in the 10" aperture is 1.003, 1.005, 
1.006 and 1.008 in BVRI. The corrected flux-ratios of 
the individual measurements to the median of the OOT 
were calculated in the manner FRcorr = 1 + S ■ (FR — 1), 
and are shown in Table |3| There is a small differ- 
ence (~2%) between the 10" and 45" flux contribution, 
thus we expect that the former measurements (0HP1.2, 
Wisel.O, TopHAT) show slightly deeper transits than the 
FLW01.2m (r) and wide-field HATNet telescopes (I). 

3. DERIVING THE PHYSICAL PARAMETERS OF THE 
SYSTEM 

We use the full analytic formula f or no nlinear limb- 
darkening given in iMandel and AgoM l|2002t) to calculate 
our transit curves. In addition to the orbital period and 
limb-darkening coefficients, these curves are a function of 
four variables, including the mass (M*) and radius (i?*) 
of the star, the radius of the planet (i?p), and the incli- 
nation of the planet's orbit relative to the observer (ip). 
Because these parameters are degenerate in th e transit 
curve, we use M, = 0.82 ± 0.031/© from IB05I to break 
the degeneracy. As regards i?*, there are two possible 
approaches: i) assume a fixed value from independent 
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TABLE 3 
The light-curve of HD 189733. 



Tel. 


Fil. 


Ntr 


HJD 


Mag 


Merr 
(mag) 


FR 


FRcorr 


FRcrr 


(mag) 


(mag) 


Qflag 


OHP1.2 


B 





2453629.3205430 


8.6062 




0.99614 


0.99612 










FLW01.2 


r 


6 


2453642.6001600 


7.1886 


0.0008 


1.02934 


1.02972 


0.00074 








HAT-5 


I 


6 


2453642.5903353 


6.7452 


0.0017 


0.99522 


0.99511 


0.00156 


-0.147 


0.0142 


G 


HAT-6 


I 


6 


2453642.6082715 


6.7357 


0.0017 


1.00397 


1.00406 


0.00156 


-0.011 


0.0160 


G 


TopHAT 


V 


6 


2453642.6042285 


7.6717 


0.0009 


0.99844 


0.99841 


0.00083 


-0.097 


0.0080 


G 


HAT-9 


I 


7 


2453644.8307479 


6.7383 


0.0019 


1.00157 


1.00160 


0.00175 


-0.015 


0.0081 


G 


Wisel.O 


B 


22 


2453678.1963150 


8.6374 


0.0013 


0.96792 


0.96779 


0.00120 









Note. — This tabic is published in its entirety (2938 lines) in the electronic edition of the paper. A portion is shown here 
regarding its form and content with a sample line for each telescope in the order they observed a transit with an OOT section. 
Column Ntr is the number of transits since the discovery data by OHP on HJD — 2453629.3. Values in the Mag (magnitude) 
column have been derived by shifting the zero-point of the particular dataset at Nt^- to bring the median of the OOT section 
to the standard magnitude value in the literature. Merr (and FRcrr) denote the formal magnitude (and flux-ratio) error 
estimates based on the photon and background noise (not available for all data). The flux-ratio FR, shows the ratio of the 
individual flux measurements to the sigma-clipped median value of the OOT at that particular Ntr transit observation. The 
merger-corrected flux-ratio FRcorr is described in detail in ^2.61 The AMext is a measure of the extinction on a relative scale 
(instrumental magnitude of reference minus image), CTmfit is the rms of the magnitude flt between the reference and the given 
frame. Both of these quantities are useful measures of the photometric conditions. Qflag is the quality flag: "G" means good, 
"C" indicates that the measurement should be used with caution, e.g. the star was marked as saturated. Fit of the transit 
parameters were performed using the HJD, FR, FRcorr and FRcrr columns. 



measurements (t l3.1(l . ii) measure the radius of the star 
directly from the transit curve, i.e. leave it to vary freely 
in the fit. Our final results are based on detailed analysis 
(fJ23 using the first approach. To fully trust the second 
approach, one would need high precision data with rela- 
tively small systematic errors. Nevertheless, in order to 
check consistency, we also performed an analysis where 
the stellar radius was left as a variable in the fit, and 
also checked the effect of systematic variations in the 
light-curves (see later in ii3.2(l . Refined cphcmerides and 
center of transit time residuals are discussed in 



3.1. The radius of HD 189733 

Because the value of the stellar radius we use in our fit 
linearly affects the size of the planet radius we obtain, we 
use several independent methods to check its value and 
uncertainty. 

First method — For our first calculation we use 
2MASS llCutri ct al. 2003) and Hipparcos photometry 
ijPerrvman et al. 1997) to find the V-band magnitude 
and V-K c olors of the st a r, and use the relation de- 
scribed in iKervella et alJ l|20n4D to find the angular 
size of the star. Because this relation was derived us- 
ing Johnson magnitudes, we first convert the 2MASS 
Ks = 5.541 ± 0.021 magnitude to the Bessell-Brett ho- 
mogenized system, which in turn is based on the SAAO 
system, a nd thus is the c losest to Johnson magnitude 
available l)CarT)enteill2005j) . We obtain a value K = 
5.59±0.05. Most of the error comes from the uncertainty 
in the J — Kg color, which is used in the conversion. 

The Johnson V-band magnitude from Hipparcos is 
7.67 ±0.01. This giv es a V-K color of 2.09 ± 0.06. From 
IKervella et al.l 1)20041) . the limb-darkened angular size of 
a dwarf star is related to its K magnitude and V-K color 
by: 



log{e) = 0.0755(1/ -K) + 0.5170 - 0.2K . 



(1) 



Given the proximity of HD 189733 (19.3 ± 0.3pc), red- 
dening can be neglected, despite its low galactic latitude. 
The relation gives an angular size of 0.36±0.02mas for the 
stellar photosphere, where the error estimate originates 



from the errors of V-K and K. The small dispersion of the 
relation was not taken i nto account in the erro r estimate, 
as it was determined bv IKervella et al ! (2004) using a fit 
to a sample of stars with known angular diameters to be 
less than 1%. Using the Hipparcos parallax we find that 
= 0.75 ± 0.05i?o. 

Second method — We also derive the radius of the star di- 
rectly from the Hipparcos parallax, V-band magnitude, 
and temperature of the star. We first convert from ap- 
parent magnitude to absolute mag nitude and apply a 
bolometric correction l|Bessell et al.l lT998). To solve for 
the radius of the star we use the relation: 



Mb = 4.74 - 2.blog 



eff,* 



(2) 



For Teff^* = 5050 ± 50 K effective temperature l|BO 
we measure a radius of 0.74 ± O.O3i?0. 

Third method - isochrones — An additional test on the 
stellar radius and its uncertainty c omes from stellar evo- 
lution models. We find from the IGirardi et al.l l)2002fl 
models that the isochrone gridpoints in the {Teg, log 5) 
plane that are closest to the observed values [Teff^^, = 
5050 ± SOif, logg = 4.53 ± 0.14) prefer slightly evolved 
models with M* « O.8OM0 and « O.79i?0. Alterna- 
tively, the Hipparcos V ~ 7.67 ± 0.01 magnitude, com- 
bined with the m~M ~ — 1.423±0.035 distance modulus 
yields My — 6.25 ± 0.04 absolute V magnitude, and the 
closest isochrone gridpoints prefer less evolved stars with 
M, « O.8OM0 and i?* « O.76i?0. The discrepancy be- 
tween the above two approaches decreases if we adopt a 
slig htly larger distance modulus. Finally, comparison to 
the lBaraffe et all (|1998l) isochrones yields « O.8OM0 
and « O.76i?0. From isochrone fitting, the error on 
the stellar radius can be as large as O.O3i?0. 



Fourth method — Recently IMasana etaLl pool cal- 
ibrated the effective temperatures, angular semi- 
diameters and bolometric corrections for F, G, K type 
stars based on V and 2MASS infrared photometry. 
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TABLE 4 

Parameters from simultaneous fit of 
transit curves. 



Parameter 



Best-Fit Value 



Rp{Rj) 

M,(Mq) 
R*{Re) 

Period (days) 
To (HJD) 



1.154 ±0.032 
85.79 ±0.24 
0.82 ±0.03'' 
0.758 ±0.016'' 
2.218573 ±0.000020 
2453629.39420 ± 0.00024 



°-From lBn,'3 

''From lMasana et all i200i 



They provide ~ among other parameters - angular semi- 
diameters and radii for a large sample of Hipparcos stars. 
For HD 189733 they derived R = 0.758 ± O.OlGi?©. 

Summary — Altogether, the various methods point to a 
stellar radius in the range of 0.74 to 0.79i?o, with mean 
value bei ng ^ Q.lQRp,. In the subsequent analysis we 
accept the lMasana et al.l l)2006j) value of O.758±O.O16iT!0. 

3.2. Fitting the transit curve 

We set the mass and radius of the star equal to 0.82 ± 
O.O3M0 and 0.758±0.016 i?©, respectively, and fit for the 
planet's radius and orbital inclination. The goodness-of- 
fit parameter is given by: 



N 

E 



(3) 



where rui is the measured value for the fiux from 
the star (with the median of the out of transit points 
normalized to one), pi is the predicted value for the fiux 
from the theoretical transit curve, and am,i is the error 
for each fiux measurement. 

For the 0HP1.2 and Wise data, where independent er- 
rors for each fiux measurement are not available, we set 
the am,i errors on all points equal to the standard devia- 
tion of the out of transit points. For the FLW01.2, HAT 
and TopHAT data, where relative errors for individual 
points are available, we set the median error equal to the 
standard deviation of the out of transit points and use 
that to normalize the relative errors. We also allow the 
locations of individual transits to vary freely in the fit. 

When calculating our transit cur ves, we use the non- 
linear limb-darkening law defined in IClare^ l|2b 00'l: 



/(r) = l-^c„(l-A."/2) 



where 



/i = cost 



(4) 



(5) 



We select the fo ur-parameter n onlinear limb-darkening 
coefficients from IClaretl l|200CI|) for a star with T = 
5000X, log{g) = 4.5, [Fe/H] = 0.0, and a turbulent 
veloc ity of 1.0 km/s. The actual parameters for the star, 
from IB051 are rather close to this: T = 5050 ± 50if , 
log{g) = 4.53 ± 0.14, and [Fe/H] = -0.03 ± 0.04. 

To determine the best-fit radius for the planet, we eval- 
uate the function over all full transits simultaneously. 





n9R l 

-0.10 -0.05 0.00 0.05 0.10 
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FLW01.2 r', N,=6 
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FLW01.2 r', binned 
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TopHAT V, N„=6 
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Fig. 1. — The five full eclipses examined in this work, with best- 
fit transit curves over-plotted. The figure in the electronic edition 
is color-coded according to the bandpass used. 

using the same values for the planetary radius and incli- 
nation. For this purpose, we employed the d ownhill sim- 
plex m inimization routine (amoeba) from iPress et alJ 
i'1992^. The full transits and the fitted curves are ex- 
hibited on Fig. ^ the transit parameters are listed in 
Tabled In order to determine the la errors, we fit for 
the inclination and the radius of the planet using the la 
values for the mass and radius of the star (assuming they 
are uncorrelated). We find that the mass of the star con- 
tributes errors of ±0.004 i?j and ±0.12°, and the radius 
of the star contributes errors of ±0.032 i?j and ±0.21°. 
Using a bootstrap Monte Carlo method, we also esti- 
mate the errors from the scatter in our data, and find 
that this scatter contributes an error of ±0.005 i?j and 
±0.03° to the final measurement. This gives us a total 
error of ±0.032 i?j for the planet radius and ±0.24° for 
the inclination. 

Our best-fit parameters gave a reduced value of 
1.23. The excess in the reduced over unity is the re- 
sult of our method for normalizing the relative errors for 
data taken at Ntr = 6, where the RMS variation in the 
data increases significantly towards the end of the data 
set, as the source moved closer to the horizon. For these 
data we define our errors as the standard deviation of 
the data before the transit, where the scatter was much 
smaller. This is justified because we know from several 
sources (night webcamera, raw photon counts) that the 
conditions were similar (photometric) before and during 
the transit, and the errors before the transit better rep- 
resent those inside the transit. This underestimates the 
errors for data after the transit, inflating the function 
accordingly. We find that when we exclude the FLW01.2 
data after the end of the transit (the FLW01.2 data con- 
tain significantly more points than any other single data 
set), the reduced for the fit decreases to 0.93. 

The results of the planet transit fit are shown in 
Tabled The value for the radius of the planet Rp = 
1.15 ± 0.03 -Rj is smaller than the B05 value (.Rp = 
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1.26 ± 0.03 and the inclination of 85.8 ± 0.2° is 

slightly larger than the lB05l value (85.3±0.1°). Altho ugh 
our errors are comparable to the errors given bv lBOSl de- 
spite the superior quality of the new data, we note that 
this is a direct result of the larger error (±0.016 instead 
of ±0.01) for the stellar radius we use in our fits. As 
discussed in H3.ll we feel that this error, which is based 
on the effective temperature and bolometric magnitude 
of the star, is a more accurate reflection of the uncer- 
tainties in the measurement of the radius of the star. We 
note that the errors are dominated by the uncertainties 
in the stellar parameters (notably i?,). 

Fitting with unconstrained stellar radius — We note that 
when we fit for the stellar radius directly from the transit 
curves (meaning we fit for the planet radius, orbital incli- 
nation, and stellar radius, but set the stellar mass equal 
to 0.82 Mq), we measure a stellar radius of 0.678 ±0.015 
Rq and planet radius of Rp = 0.999 ± 0.026i?j. The er- 
rors for these measurements are from a bootstrap Monte 
Carlo analysis, and represent the uncertainties in our 
data alone. To obtain the formal errors, we incorpo- 
rate the error from the mass of the star and find er- 
rors of ±0.017i?Q and 0.029i?j, respectively. This means 
that our data prefer a significantly smaller stellar radius 
(and a correspondingly smaller planet radius) than our 
estimates based on temperature, bolometric magnitude, 
and V-K colors alone would lead us to expect, or a radius 
smaller than the O.82Af0 stellar mass implies. 

With many more points (869 as compared to ~ 100 in 
the other data-sets) and lower photon-noise uncertain- 
ties, the FLW01.2 data dominate the fit of Eq. How- 
ever, we repeated our fit with and without these data, 
and found that the best-fit radius for the star decreased 
only slightly (to 0.666 Rq) when the FLW01.2 data were 
excluded from the fit. Thus, our I, V, and B-band data 
independently yield values for the stellar radius similar 
to those implied by the FLW01.2 data. 

The effect of systematic errors — The minimization for- 
mula (Eq. O assumes independent noise, and the pres- 
ence of covariance in the data (due to systematics in the 
photometry) means that too much weight may be given 
to a dataset having small formal errors and a great num- 
ber of datapoints (e.g. the FLW01.2 data) compared to 
the other independent data sets (e.g. other telescopes and 
filters). This is especially a concern when the datasets 
yield different transit parameters, and one needs to es- 
tablish whether this difference is significant. In order to 
follow-up this issue, we repeated the global fit by assum- 
ing that the photometric systematics were dominant in 
the error budget on the parameters - as suggested by our 
experience with milli-magnitude rapid time-series pho- 
tometry. We estimated the amplitude of the covariance 
from the variance of 20-minute sliding averages on the 
residuals around the best-fitting transit light-curve for 
each night following the method of Pont (2005). The fit 
was repeated using these new weights (listed in Tabled 
(Jsys), and the resulting parameters (planetary radius, in- 
clination) were within 1% of the values found assuming 
independent noise. The dispersion of these parameters 
from the individual nights were found to be compatible 
with the uncertainties due to the systematics. The am- 
plitude of the systematics is also sufficient to account for 
the difference in the best-fit stellar radius if it is left as 
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Fig. 2. — The ten partial eclipses examined in this work, with 
best-fit transit curves over-plotted. The eclipses are listed sequen- 
tially by date, from top left to bottom right. These eclipses were 
not used in the fit for the planet radius, inclination, stellar mass, 
and stellar radius. The figure in the electronic edition is color- 
coded according to the bandpass used. 

a free parameter. Therefore, with the amplitude of the 
covariance in the photometry determined from the data 
itself, we find that the indications of discrepancy between 
the different data sets and with the assumed primary ra- 
dius are not compelling at this point. 

3.3. Ephemerides 

The transit curves derived from the full-transits for 
each bandpass were used in turn to calculate the 
ephemerides of HD 189733 using all transits that have 
significant OOT and in-transit sections present (for ref- 
erence, see Table |21l. For each transit (full and partial), 
the center of transit Tc was determined by minimiza- 
tion. Partial transits with the fitted curve overlaid are 
exhibited on Fig.|2| Errors were assigned to the Tc values 
by perturbing Tc so that increases by unity. The indi- 
vidual Tc transit locations and their respective errors are 
listed in Table [S] The typical timing errors were formally 
of the order of 1 minute. This, however, does not take 
into account systematics in the shape of the light-curves. 
The errors in Tc can be estimated from the simultane- 
ous transit observations, for example the Ntr = 6 event 
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Telescope Nt 



Tc 

(HJD, days) 



<^HJD 

(days) 



(O-C) 
(days) 



(O-C) 
■^HJD 



OHP1.2 

OHP1.2 

OHP1.2 

FLW01.2 

HAT-5 

HAT-6 

TopHAT 

HAT-9 

HAT-9 

HAT-5 

TopHAT 

HAT-5 

Wise 

TopHAT 

HAT-9 




4 
5 
6 
6 
6 
6 
7 
16 
19 
19 
20 
22 
24 
29 



2453629.39073 
2453638.26885 
2453640.48706 
2453642.70592 
2453642.70641 
2453642.70649 
2453642.70536 
2453644.92720 
2453664.89287 
2453671.54999 
2453671.54849 
2453673.76725 
2453678.20080 
2453682.63715 
2453693.73327 



±0.00059 
±0.00067 
±0.00174 
±0.00022 
±0.00092 
±0.00049 
±0.00048 
±0.00111 
±0.00108 
±0.00113 
±0.00096 
±0.00072 
±0.00050 
±0.00100 
±0.00090 



-0.0035 
0.00035 
-0.0000079 
0.00029 
0.00077 
0.00085 
-0.00028 
0.0030 
0.0015 
0.0029 
0.0014 
0.0016 
-0.0020 
-0.0028 
0.00045 



-5.9 
0.53 
-0.0045 
1.24 
0.84 

1.7 
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Fig. 3. — The residuals calculated using the period and Tq 
derived in this work. The dashed lines are calculated from the 
uncertainties in the measurements of P and To . The figure in the 
electronic edition is color-coded according to the bandpass used. 



and partial eclipses examined in this work 
transits and O-C residuals for each eclipse. 



was observed by the FLW01.2m, HAT-5, HAT-6 and 
TopHAT telescopes (Table and the rms of Tq around 
the median is ~ 50 seconds, which is in harmony from 
the above independent estimate of 1 minute. We ap- 
plied an error weighted least square minimization on the 
Tc ~ P-Ntr+E equation, where the free parameters were 
the period P and epoch E. The refined ephemeris values 
are listed in Table 0] They arc co nsistent with both those 
derived by B05 and by Hcbrard fc Lecavelier Pes Etangj 
l)2006fl using Hipparcos and 0HP1.2 data, to within lo- 
using our error bars. 

We also examined the Observed minus Calculated (O- 
C) residuals, as their deviation can potentially reveal the 
presence of moons or additional planetary companions 
l|Holman fc Murravl l20(i5l: lAgol et al.ll2nf)5(l . The O-C 
values are listed in Table 0] a nd plotted in Fig. El Usin g 
the approximate formula fromlHolm an fc Murravl l)2005f) . 
as an example, a 0.15A/j perturbing planet on a circular 
orbit, at 2 times the distance of HD 189733b (P « 6.3'^) 
would cause variations in the transit timings of 2.5 min- 
utes. The radial velocity semi-amplitude of HD 189733 
as induced by this hypothetical planet would be 19m/s, 
which would be barely noticeable (at the la level) from 
the discovery data having residuals of 15m/s and span- 
ning only 30 days. 

A few, seemingly significant outlier points on the O- 
C diagram are visible, but we believe that it would be 
premature to draw any conclusions, because: i) the error- 
bars do not reflect the effect of systematics, and for ex- 
ample the Tc of the Ntr = OHP discovery data moved 
by ^ 5 minutes after re-calibration of that dataset, ii) 
all negative O-C outliers are B or V-band data, which 
is suggestive of an effect of remaining color-dependent 
systematics. The significance of a few outliers is further 
diminished by the short dataset we have; no periodicity 
can be claimed by observing 2 full and 9 partial transits 
altogether. 

According to the theory, the nature of perturbations 
would be such that they appear as occasional, large out- 
liers. Thus, the detection of potential perturbations also 
benefits from the study of numerous sequential transits. 



for example, the MOST mission with continuous cover- 
age a nd uniform da ta would be suitable for such study 
We also give the number of elapsed! WalkereLalJl2003(). draw the attention to the 

importance of observing full transits, as they improve 
the Tc center of transit by a significant factor, partly 
because of the presence of ingress and egress, and also 
due to a better treatment of the systematics. 

If a planet is perturbed by another planet, the transit- 
time variations At arc proportional to the period P of the 
perturbed planet (Holman & Murray 2005). Although 
HD 189733b is a relatively short period (2.21'^) planet 
compared to e.g. HD 209458b (3.5'*), it is a promising 
target for detecting transit perturbations in the future, 
because the mass of the host star is low and At oc 1 /AI^, , 
plus the deep transit of the bright source will result in 
very precise timing measurements. Observations span- 
ning several months to many years may be needed to 
say anything definite about the presence or absence of a 
periodic perturbation. 



4. CONCLUSIONS 

Our final values for the planet radius and orbital in- 
clination were derived by fixing the stellar radius and 
mass to indep e ndent ly determined values from lB05l and 
iMasana et al.l l)2006j) . We analyzed the dataset in two 
ways: by minimization assuming independent errors, 
and also by assuming that photometric systematics were 
dominant in the error budget. Both methods yielded the 
same transit parameters within 1%: if we assume there 
is no additional unresolved close-in stellar companion to 
HD 189733 to make the transits shallower, then we find 
a planet radius of 1.154 ± 0.032i?j and an orbital incli- 
nation of 85.79 ±0.24° (TableEj). The uncertainty in Rp 
is primarily due to the uncertainty in the value of the 
stellar radius. 

We note that the TopHAT V-band full and partial 
transit data, as well as the Wise partial B-band data 
appear slightly deeper than the best fit to the analytic 
model. The precision of the dataset is not adequate to 
determine if this potential discrepancy is caused by a real 
physical effect (such as a second stellar companion) or to 
draw further conclusions. 

When compared to the discovery data, the radius de- 
creased by 10%, and HD 189733b is in the mass and ra- 
dius range of "normal" exoplanets (Fig. The revised 
radius estimate is consistent with structural models of 
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a higher density of Igcm as compared to the former 
~ 0.75gcm^'^. The smaller planetary radius increases 
the le^im brightness temperature T^^'^^) = 1117 ± 42K 
of Demine; et al. (2006) to 1279 ± 90K, which is slightly 
larger than that of TrES-1 and HD 209458b. 

We also derived new ephemerides, and investigated the 
outlier points in the 0-C diagram. We have not found 
any compelling evidence for outliers that could be due to 
perturbations from a second planet in the system. We 
note however, that due to the proximity and brightness 
of the parent star, as well as the deep transit, the system 
is well suited for follow-on observations. 



Fig. 4. — Mass-radius diagram depicting known transiting ex- 
oplanetSj plus Saturn and Jupiter (for comparison). Both the dis- 
covery : lifl.";. small filled square) and revised radius (large filled 
square) for HD 189733 is shown. Sources for the mass and ra- 
dius values (in this order) are l isted a fter the name of the planet; 
HD 209458b; ILaughlin eFal? f2003), Knutson et. a l. 2006, in 
prepar ation, HD 189733: BOS; this pa per, OGLE lllb:IPont et al.l 
120041) . OGLE lOb: Konacki et al. 1200^. Igantos et alj 11200611. 
OGLE 132b: Moutou et al. (2004), OGLE 56b: ITorres etaL 

fl04), OGLE 113 b: Bouchv et al.. L2n04l TrES-1: ILaughlin et al. 
105) HD 149026: ISato et alj 120051) . ICharbonneau et alJ ll2006l " 

hot Jupiters that include the effects of stellar insolation, 
and hence it does not require the presence of an addi- 
tional energy source, as is the case for HD 209458b. On 
the mass-radius diagram HD 209458b remains an out- 
lier with anomalously low density. We note that t he pa- 
rame t ers of OGLE-lOb ar e still debated (Konack i et al.l 
120051 iHolman et all 120061: iSantos et aL 2006), but ac- 
cording to the recent analysis of ISantos et .aL. {2DQS)i 
it also has anomalously low density. With its revised 
para meters, HD 189733b is quite similar to OGLE-TR- 
132b ijMoutou et alJl2004D . The smaller radius leads to 
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